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Abstract 

The Heisenberg model, a quantum mechanical analogue of the Ising model, has a large ground 
state degeneracy, due to the symmetry generated by the total spin. This symmetry is also 
responsible for degeneracies in the rest of the spectrum. We discuss the global structure of the 
spectrum of Heisenberg models with arbitrary couplings, using group theoretical methods. The 
Hilbert space breaks up in blocks characterized by the quantum numbers of the total spin, S 
and M, and each block is shown to constitute the representation space of an explicitly given 
irreducible representation of the symmetric group Sn, consisting of permutations of the N spins 
in the system. 

In the second part of the paper we consider, as a concrete application, the model where each 
spin is coupled to all the other spins with equal strength. Its partition function is written as 
a single integral, elucidating its A^-dependence. This provides a useful framework for studying 
finite size effects. We give explicit results for the heat capacity, revealing interesting behavior 
just around the phase transition. 
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1 Introduction 



The Ising model is one of the best known models of interacting spins. In its simplest form, 
with nearest-neighbor interactions only, the Hamiltonian is 

Hlsing J ^ ^ "^jSj ) (1) 

(ij) 

where the spins take the values +1 (up) or —1 (down) and the summation is over nearest- 
neighbor pairs (ij). The ground state of the model depends on the sign of the coupling pa- 
rameter J. If J > (ferromagnetic behavior) the spins tend to align themselves with their 
neighbors, and the ground state is the configuration with either all the spins up or all down. 
For J < (antiferromagnet) the ground state is the checkerboard-like Neel state, with all the 
spins up at the even site^ and down at the odd sites or the other way round. The twofold de- 
generacy of the ground state in both cases reflects the spontaneous breaking of the Z2 symmetry 
of flipping all the spins. 

The Heisenberg model is a quantum mechanical analogue of the Ising model (for an intro- 
duction see e.g. ref. m, an early discussion can be found in ref. P]). The spectrum of this model 
is not as well understood, in particular the antiferromagnetic ground state is not known. In the 
Heisenberg model, however, there is a global rotation symmetry, associated to the operator of 
the total spin, which is very helpful in determining the global structure of the spectrum. This 
applies to the model in its most general form, i.e. with arbitrary coupling between each pair of 
spins. 

In the first part of this paper we discuss the global spectral structure of Heisenberg mod- 
els from a group theoretical point of view. We start with a pedagogical introduction to the 
Heisenberg model (section (^), demonstrating the role of the total-spin operator and writing 
down explicitly the degenerate ferromagnetic ground states. In section |^ we consider the con- 
sequences of the global rotational symmetry for the rest of the spectrum. It is shown that 
the spectrum breaks up into blocks labelled by the quantum numbers of the total spin and 
characterized by an explicitly given representation of Sn, the symmetric group on N elements. 
Parts of this presentation are well known, but our approach provides an attractive alternative 
point of view, which exposes clearly the ^Ar-representation structure of the spectrum. 

In section ^ we consider the isotropic infinite range Heisenberg model, with all the spins 
equally coupled. In this model, knowledge of the global structure of the spectrum is sufficient 
to write the partition function in the form of a single integral which shows a simple dependence 
on A^. This result is particularly interesting for studying finite size deviations from the model's 
thermodynamic (A^ 00) behavior. As a concrete application we present calculations for the 
heat capacity. 



site i = (ii, Z2, ■ • ■ , id) is called even (odd) when ii + 12 + ■ ■ • + irf is even (odd). 
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2 The Heisenberg model 

In the Heisenberg model, the spins are operators in a Hilbert space. The Hamiltonian is 

H = -JY.s^-s, (2) 

{ij) 

where the components s^,s^, sf of a single spin Si constitute a set of generators of the rotation 
group SU(2). The spins satisfy the commutation relations 

[s^, 4] = t6,,e^'^sl (3) 

For a single spin, the states are labelled by the quantum numbers {s,m), determined by the 
eigenvalues s{s + 1) and m of and s^, respectively. We shall restrict ourselves to the spin-i 
case, which bears a close resemblance to the Ising model. We take = where the a"" are 
the Pauli matrices 

in terms of which the up state is ( q ), with quantum numbers (s = |,m = |), and the down 

state is ( "i* ) with (s = = — i). 

It is convenient to write the Hamiltonian as 

H = -JY^^t^j + + s-sl), (5) 

in terms of the usual ladder operators 



= ± is^ (6) 







with the properties 



(7) 



The s^s| term in (|^) has the same effect as the Ising Hamiltonian (|1|). It tends to align 
the neighboring spins (anti)parallel in the (anti)ferromagnetic case. The term between 

brackets in (|^), sometimes called the fluctuation term, is special to the quantum case. It gives 
zero when acting on parallel neighboring spins, so that the ferromagnetic Ising ground states 
with all spins up or all down are also ground states of the Heisenberg ferromagnet. Its effect on 
a pair of opposite spins, however, is to flip both the spins, so that the Neel states are not even 
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eigenstates of the Heisenberg antiferromagnet. For a two-spin system the solution is of course 
to take the antisymmetric (singlet) combination 



ir» 



but for a system of many spins the antiferromagnetic ground state is not known. 

However, there is more to the ferromagnetic Heisenberg model as well. Its ground state 
degeneracy is not just twofold but (A^ + l)-fold, where is the number of spins in the system 
(see 0, ^). The reason for this large degeneracy is that the total spin operator S commutes with 
the Hamiltonian, providing us with the quantum numbers S and M. Since the components 
gx,y,z mutually commute, one can create states with M = —S + 1, —5* + 2, . . . , S, 

degenerate in energy, by repeated application of 5*+ on a state with M = —S. (Of course, the 
degeneracy is lifted when the global symmetry is broken by an external magnetic field, through 
the addition of an interaction term 



H 



B 



(9) 



to the Hamiltonian. This causes a 'fanning out' of the energy levels, the familiar Zeeman effect.) 

It is amusing to see what the non-trivial ground states actually look like, using this ladder 
procedure. Starting from the state with all spins down one obtains the states 



|0) 
|1) 

|2) 




i) + un---i) + ■■■ + u---iT)) 



TTii---i) + ITin---i) + ■■■ + U---iiTT)) 



|N) 



TT---T). 



The state |K) is the properly normalized symmetric sum over the 
up and N — K spins down. Its quantum numbers are 



(10) 

states with K spins 



S = -N, 
2 ' 



M 



K--N , 
2 ' 



in terms of which 



S^\K) 



{S ± M + 1){S T M) |K±1) 



(12) 



At first sight it may seem surprising that the states |K), containing antiparallel spin pairs, 
are ground states of the ferromagnet. A particularly interesting one is |^N), a linear combi- 
nation of states with as many up as down-spins (even the Neel states contribute to it!). More 
generally, for finite each of the 2^ basis vectors |±, it, ■ ■ ■ , ±) (where + =1, — =|) has a 
nonzero projection on one of the (A^-|- 1) ground states of the ferromagnet. This indicates that 
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this set of basis vectors is not very suitable to describe the spectrum of the Heisenberg model. 
The reason is of course that the up and down states are defined with respect to a pre-defined 
^-axis, obscuring all rotational symmetry. 

To illustrate the relation with global rotational symmetry, consider for example the expec- 
tation value of the two point function {s^ s^) in the different ground states. For the states |0) 
and |N) one finds 

(^■^■) = 7' (13) 



but the other ground states give 



{K|s-sJ|K> = i('l-4 



^'/"-^'(l-^,)^ . (14) 



Summation over i and j gives (|A^ — K)'^ = M^, confirming the eigenvalue of |K). The 



average of eq. (|T^) over the ground states |K) is 



av(s-S-)av=^(l + 25.,)- (15) 

For i j, this is one third of the value | it attains in the states |0) and |N), as expected from 
rotational invar iance. 

The rotational symmetry generated by the total spin is a general feature of Heisenberg 
models, it is not typical of the nearest-neighbor model with equal interaction strengths we 
have been considering so far. In fact, the entire spectrum breaks up in blocks labelled by the 
quantum numbers of the total spin, each block being characterized by a certain irreducible 
representation of Sn, the symmetric group, or permutation group, on elements. This is the 
topic of the next section. 

3 Structure of the spectrum 

We will start out with a general spin-| Hamiltonian for a system of spins, 

H = -J2 Jij ■ (16) 
1^/1 \ 

= -2 E + 2^^ts- + s-sl)] . (17) 

The arbitrary interaction strengths Jjj depend on i and j now and are not restricted to nearest- 
neighbor pairs. We assume 

Jij = Jji , Jii = foralH,j (18) 

without loss of generality. 
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As a consequence of the general form of this Hamiltonian, neither the dimensionahty of the 
problem nor the boundary conditions are specified a priori. All the information is contained in 
the set of coupling parameters Jij. In the nearest-neighbor model, for example, the dimension is 
determined by the number of nearest neighbors of a given spin, which is the number of non-zero 
couplings Jij for fixed i. Similarly, (anti)periodic or free boundary conditions can be imposed 
by choosing the Jij appropriately. 

It is not difficult to show that H is invariant under global rotations. The A^-spin Hilbert 
space is the tensor product of A^ single-spin Hilbert spaces. In terms of this tensor product, 
the operator for the total spin has the form 

S = y^Si = s®/®---®/ + /®s®/®---®/ 

i (19) 
+ ■■■ + J®/® ■■■®/®s , 

where the i*'' factor in each term acts on the i^^ spin. The commutator of a single-spin compo- 
nent with the Hamiltonian is 

j 

From this it is clear that 

[S,H] = (21) 

since each term on the right hand side of eq. ( pOD is canceled by a similar contribution from 
[s'^,H]. This means conservation of total spin 5*. Its components satisfy the commutation 
relations 

[^«, S''] = ie"''" S" (22) 

so the states in the spectrum carry the quantum numbers S and M. 

We will follow a different approach, though. The Hamiltonian can be interpreted as an 
element of the group algebra of S^, the symmetric group on A^ elements. From this point of 
view, the Hilbert space is the representation space for an 5'7v-representation on which H acts. 
This space will be shown to break up into blocks and S and M will arise as natural labels on 
these blocks. The states in each block constitute an irreducible representation of Sn which is 
given explicitly. 

First, we will show that the Hamiltonian can be viewed as an element of the group algebra 
of Sn- The group algebra of a group G is the set of formal sums J^gf^c^g ' 9i where the Xg 
are numbers, with the obvious addition and multiplication rules. The group algebra plays an 
important role in the representation theory of finite groups. Now consider a term {ij) in if, 

^fo-> = sis] + i(s+s7 + sis+) (23) 

(ignoring the prefactor for the moment). Its action only depends on the spins i and j of the 
state in Hilbert space on which it acts. Denoting the values of these spins by |^, j), we have 

^fo->ITT) = JiTT), 
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ii) = 




Ti) = 




it) = 



^1 ii) 5 

i(-l|Ti) + 2| it)), 

i(-lUT) + 2|Ti)). (24) 



This can be written in a general form as 



H{^J)\^.j) = \{-l\i,j) + 2\j,i)) 

= 1(2(2j)-/)Kj), (25) 

where {ij) G Sn is the 2-cycle permutation or transposition interchanging the spins at sites i 
and j. The contribution of the identity operator /, which is the unit element of Sn, can be 
removed by adding a constant to the Hamiltonian (|l^), and will be neglected. The Hamiltonian 
is a linear combination of terms of the form (^), with coefficients proportional to Jjj, 

H = -WJ^^i^3)■ (26) 

^ {ij) 

As such, it is an element of the group algebra of S^- Note that it is restricted to the subspace 
of this algebra spanned by the conjugacy class of transpositions (2-cycles). In the calculation 
of the partition function, however, one needs the exponentiated Hamiltonian exp[— which 
is not restricted to this conjugacy class. 

The 2^-dimensional Hilbert space of the spin model can thus be viewed as the representation 
space of a representation of Sn- It is clear that this representation is reducible. Any permutation 
leaves the numbers of up and down-spins unchanged, so the representation space decomposes 
into a sum of invariant subspaces Wk of dimension ( x )' corresponding to the sectors of 

different S^ eigenvalues M = K — \N . Each basis element of the subspace Wk can be labelled 
by the K sites where the up-spins are located, 

= {() = |0)} 
W, = {(l),(2),(3),...,(Ar)} 

Wk = {iii,i2,---,iK)\l<ii<i2<---<iK<N} 



Wn = {(l,2,3,...,Ar) = |N)}. (27) 
By taking the symmetric sum of all the states in Wk, for a fixed value of K, one recovers the 



state |K) of eq. (pl)|), which is a singlet under Sn- 

Indeed, the representation subspaces Wk are not irreducible under Sn either. Not surpris- 
ingly, they decompose into irreducible sectors with different total spin quantum number S. This 
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decomposition, described nicely in Wigner's book 0, can be summarized as follows. Within 
each subspace Wk (for K < ^N) there is a subset of states which transform as Wk-i- This 
gives a decomposition of Wk as 

Wk = Wk-1 © Dk (28) 

where Dk can be shown to be an irreducible subspace. Subsequently, one applies the same 
decomposition for Wk~i and so on until one ends up with the full decomposition 

Wk = Do®Di®D2®---®Dk. (29) 

Finally, the decomposition of Wk for K > is found by noting that 

Wn-k ^ Wk . (30) 

Before discussing the representation spaces Dk, let us say something more about this re- 
cursive decomposition of Wk- Since Wk-i can be regarded as a subspace of Wk according to 
eq. (|28|), it is interesting to see how the corresponding states in the two spaces are mapped on 
each other. For every state with K — 1 up-spins we want to determine the corresponding state 
in the space spanned by basis vectors with K up-spins. Let IiIj^k-i)) = (h, "^2, ■ ■ ■ , ix-i) be an 
arbitrary state in Wk-i- Then the corresponding state \'iIj{k)) G Wk is given by the sum over 
the — -|- 1 states (ji, j2, • • • ,3k) G Wk satisfying the condition that the set of numbers 
{ia, a = 1, . . . , K} consists of the numbers {z/3, /? = 1, . . . , K — 1} supplemented by one other 
number. In other words, \'4'{k)) is the state given by the sum over all the states for which all 
the spins ia as well as one additional spin are up-spins. This is exactly the state obtained by 
acting on \ip{^K-i)) with the ladder operator 5*+. 

In terms of Young tableaux, the decomposition of Wk can be visualized as follows (see 
ref. Q). The reducible representation of Sn on the space Wk is a permutation representation 
corresponding to the so-called skew Young diagram 

N-K K 



which decomposes into the following Young diagrams for irreducible representations 



N N-l N-2 

, ^ . , ^ ^ , 

1 1 1 1 1 1 e u ' ® 



N-K 



e ■■• e 



K 



(32) 
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The irreducible representation^] Dk is identified as 



N-K 



D 



K 



(33) 



K 



Schematically, the Hilbert space decomposes as follows under S_ 



N'- 



Do 










^0 






Do 


Di 








^0 


Di 




Do 


Di 


D2 






Do 


Di 


D2 


Do 


Di 


D2 ■ 


■ ^i(7V-l) 


and 


Do 


Di 


D2 


Do 


Di 


D2 ■ 


■ ^l(^-l) 




















Do 


Di 


D2 


Do 


Di 


D2 






Do 


Di 




Do 


Di 








Do 






Do 

















Di 



N 



(34) 



for odd and even A^, respectively. In this diagram, the quantum number S increases from right 
to left, while M runs vertically. The dimension of Dk is ^ ^ — -1^ ^, with the convention 

^ = 0. It is easily checked that the dimensions of the subspaces in this decomposition add 
up to 2^. In the limit — > 00, Dij^ has dimension 2^ ■ A/N^/2'kN, while the largest is 

the one for K = — |"\/iV, whose dimension is \J N/e times as large. 

This ends our discussion of the global spectral structure. To investigate the spectrum in 
more detail one has to consider the representation theory of the S^q group algebra element H 
of eq. ( p6D in each of the representations D^- This will depend on the specific form of the 
Hamiltonian, i.e. on the couplings Jij. In practical applications, these couplings will not be 
arbitrary as in eq. ([T6| ) but related by symmetries. For example, the nearest-neighbor model 
on an L'^ lattice with periodic boundary conditions and equal couplings Jij = J for all nearest- 
neighbor pairs, has additional translational, rotational and reflection invariance properties. The 
exploitation of these extra symmetries, regarded as a subgroup of Sn = S^d, will be important 
for solving the model. 



^We use the same notation Dk for the representation itself and for its representation space. 
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4 The isotropic infinite range Heisenberg model 



As an example, we shall consider the model in which each spin interacts with every other spin 
with equal strength^. This model can be viewed as a system of N spins on the vertices of an 
(A^ — l)-dimensional simplex (the higher-dimensional generalization of a tetraeder), interacting 
along its edges. With this picture in mind, the model is effectively infinite-dimensional for 
N oo. The ferromagnetic version of this model was studied by Kittel and Shore ^ a long 
time ago. They derived expressions for the N oo limit and presented numerical results for 
various finite N. A remarkable finding was that the phase transition develops very slowly as 
N is increased. 

Here, instead of taking the N oo limit at the beginning, we first derive an expression 
for the partition function in terms of a single integral, elucidating its A^-dependence. Large- A^ 
results can be recovered subsequently by calculating this integral in the saddle-point approxi- 
mation. We will briefly mention the anti-ferromagnetic case too. 

The partition function is 

^ Tr exp[-(3H] 



with Hamiltonian 
1 

~ 2 



N 



H 



-js- 



Tr/ 



+ ^4iV(^iV + l) 



(J>0) 



(35) 



(36) 



where a constant has been added to normalize the ground state energy to zero, and 13 = l/ksT . 

For this Hamiltonian, the energy of a state is entirely determined by its quantum num- 
ber S. Calculation of the partition function thus reduces to counting the number of states 
with this quantum number. The states in the representations Dk in the scheme (^4]) have 
S = — K, and there are 25" + 1 such blocks, each consisting of 



2 



N 
K 



N 
K - 1 



N 

-{S- 



1) 



states. Hence one gets 



exp 



-PJN{N + 2) 
X E(25 + l)< 



1 

2^ 



N 



IN-S 




exp 



-(3JS{S+1) 



(37) 



adopting the convention 

for even A^ and over half integers S = |, |, 
Z can be rewritten in terms of an integral 



0. Here the summation runs over integers S = 0,1, . . . , 

, for odd A^. In the appendix it is shown that 



2 exp 



exp 



--a(Ar + l) 



d_ 

da 



drj 



^^2n/a{N + l 
{N + l)dr] 



exp 



1 2 



cosh ar] 



N+l 



^2n/a{N + 1) 



Tjianhai] exp 



cosh ar] 



N+l 



^For a discussion of the Ising analogue of this model see e.g. ref. 0. 



(3^ 



(39) 
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where 

a = ^(3J{N + l). (40) 

The advantage of these expressions is that they exhibit clearly how Z depends on N. It is 
interesting to note that the integral occurring in eq. (pH]) is very similar to the analogous 
expression for the corresponding model for classical spins, see e.g. ref. 0. 

For large A^, t] can be interpreted as the order parameter defined by the spin fraction, in 
the sense that 

2, / S{S+1 



(up to a correction). Here the expectation value on the left hand side is calculated with 
eq. ( P^D and the right hand side with eq. {^7\). This means that our rj is essentially the same 
as the one defined by Kittle and Shore 0, even though their rj is a discrete quantity taking 
values in the interval [0, 1]. In fact, from eq. (|39| ) we recover their result that for a = 1 (i.e. 
T = Tc) the most probable value of rj is (6/A^)^'^, and below we shall see ( ^31) that a similar 
correspondence holds in the general case {a ^ 1). 

We shall focus on the large- limit of the partition function and the heat capacity Cy- For 
comparison we have displayed Cy for N = 2000 in fig. 2, which is adapted from ref. [§. Guided 
by the prominent appearance of the number (A^ + 1) in eq. ( p8D we perform an expansion 
in 1/(A^ + 1), keeping a = \(3J{N + 1) fixed; the corresponding expansion in at fixed 

ao = \[3JN , is easily derived from this. For (A^ + 1) large, the integrand in eq. ( pSf ) is dominated 
by the value of r] at which the exponent 

— -7]^ + log cosh at] (42) 

is maximal. This value fj is given by the familiar condition 

tanh ar] = f/ , (43) 



signalling a phase transition at a = 1. Recall ( ^OD that a is inversely proportional to the 
temperature, so a = T^/T. 

If a < 1 the only solution to (^) is f/ = 0. For a > 1, two additional solutions appear, and 
the integrand is maximal for these non-trivial values of i), see fig. 1. A Taylor expansion of the 
integrand around of eq. ( pSj) around rj = rj now gives rise to an expansion in 1/(A^ + 1). We 
shall give results to lowest non-trivial order in 1/A^ only, although higher order corrections are 
obtained easily. For a < 1 one finds 



2 exp 



exp 



ia(Ar+l) 



#^fl + ^f-)) (44) 
da v^r^ V \NJJ ^ ' 

' '^^Oa)) , (45) 



(l-a)3/2V VA^. 

from which it follows that the lowest order correction to the value zero of the heat capacity per 
spin 
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in the high temperature phase is (for small r 



l-a^ 



For a > 1, expansion of the integrand around ±77 gives 



3 1 



1 



-a{N + l) 



Z = A exp 
which can be evaluated using 



d (exp 
da 



dfj 

da 



1 ~2 



cosh a?7 j 



N+l 



1 - a(l - 772) 



For small Irl 



1 - 0(1-772) 

^— i this gives for the heat capacity in the low temperature phase 



(47) 



(48) 



(49) 



(50) 



50D seems compatible with the finite size data of ref. P], 

|r| > 



This behavior of the heat capacity 
see also fig. 2. 

These results are valid for a not too close to 1, in fact we must have |a — 1| 
otherwise the derivations of eqs. (|45|) and (|48|) and the subsequent results for the heat capacity 
lose validity. In the thermodynamic limit N ^ 00, this means that only the phase transition 
point a = 1 itself is not described by these formulas, but for finite there is a finite interval 
of a- values around a = 1 in which the theory behaves qualitatively different. The width of this 
interval is O and can be regarded as a measure of the "width of the phase transition" for 

finite N (for convenience we keep using the term 'phase transition' although strictly speaking 
there is no phase transition in a finite system). 

This region very close to the phase transition appears to show interesting behavior. Consider 
expression (|39| ) for the partition function with r = = , where < |5| -C 1. The leading 
terms in the exponent of the integrand are 

6 



(iV + 1) 



=77 



12 



4 4 
a r] 



(51) 



which both appear to be of leading order in A^, after making the substitution 077 — > y/{N+iy^^. 
Since I (5 1 -C 1, however, we can consider the 5-dependent exponential as a perturbation. We 



find 
Z = 



exp 



exp 



X 



X 



--a{N + l) 
-ia(iV + l) 

6(i2)-i/^r(- 



1 + 



(AT +1)3/4 

ay/2na 

(AT +1)3/4 



dyy'^ exp 



12' 



O 



1 



(52) 



ay/^Tca 

i<5 6(12) V4r(^) + 6(i2)3/^r(I) - 6(12)^/^r(^) + 0(5') 



(53) 
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This leads to 




2M + (^(0 



0.61 



+ higher order corrections 



1(54) 
(55) 



Thus we have exphcitly found the value of the heat capacity and its slope at the phase transition 
in a finite system. The fact that the slope is proportional to ^/N confirms the earlier remark 
that the width of the phase transition is of order The behavior (|55D appears to be in 

qualitative and quantitative accordance with the numerical results for Cy of Kittle and Shore 
[§], see fig. 2 (recall — r ^ a — 1). 

For completeness we mention the antiferromagnetic case (see eq. (^) in the appendix) as 
well, although it is not very interesting. It might seem appropriate to redefine the normalization 



of the Hamiltonian (36) such that the antiferromagnetic ground state energy becomes zero, but 
we stick to the original normalization for continuity reasons. The integrand is maximal at f/ = 
and the partition function becomes (here a = —jPJ{N + 1) > 0) 



exp 



-a{N + l) 



[l + a)3/2 



(56) 



which is just the ferromagnetic result (|45|) with a —a. The partition function has continuous 
behavior around a = so we can regard the antiferromagnetic case as an extension of the high 
temperature phase of the ferromagnet. 



5 Conclusion 

In the first part of this paper we have discussed Heisenberg models with arbitrary couplings. 
The global structure of the spectrum of such models is determined by the well known global 
rotational symmetry provided by the total spin, which causes the spectrum to break up into 
blocks of states labelled by the quantum numbers 5* and M of the total spin. The group theo- 
retical approach pursued here shows that the states within one such block transform according 
to explicitly given irreducible representations of the symmetric group Sn- 

In an attempt to diagonalize the Hamiltonian in the representations D^, one might try to 
proceed by applying representation theory for the group algebra, in terms of the Schur func- 
tions for example. In practice, one is usually interested in Heisenberg models satisfying extra 
relations between the coupling parameters Jij, such as the nearest-neighbor model with equal 
couplings between each pair of neighboring spins. Consideration of the additional symmetries 
following from these relations, regarded as subgroups of Sjy, is probably essential. Another 
possible extension of this work would be a generalization to spins in higher-dimensional repre- 
sentations of SU(2). There, one expects the appearance of S^- representations other than the 
ones discussed here. Furthermore, it would be interesting to give a description of spin waves in 
the present framework. 

In the second part of this paper we studied the infinite range Heisenberg model with equal 
couplings between all the N spins. We have presented a compact formulation of the partition 
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function which clearly exhibits the dependence on and allows for a straightforward study of 
deviations from the N ^ oo behavior. As an example, we calculated finite size corrections to 
the heat capacity. These calculations reveal interesting behavior close to the phase transition 
and are in quantitative agreement with numerical data for finite systems. 



A Rewriting the partition function 



In this appendix we present the derivation of the integral representation (| 
function eq. 



of the partition 



Z = exp 

X 5^(25 + 1 



pJN{N + 2] 



1 

2^ 



N 



N 



lN-{S + l) 



exp 



-/3JSiS+l) 



(57) 



We will concentrate on the ferromagnetic case, J > 0. After that we will summarize the results 
for the antiferromagnet and briefly discuss the inclusion of an external field. 
The negative terms in ( |57D can be rewritten using the equality 



-(25 + i; 



A^ 

In- (8 + 1) 

= {2S' + 1) ( 



exp 



-f3JS{S + l) 



N 



exp 



-f3JS'{S' + l) 



{5t 



where S' = —{S + 1) and we have used that 
the rather elegant expression 



N 

\n + S' 



-N^-~ s' ) • ^'i- (113) ^ow reduces to 



Z = exp 



'-PJ{N+1] 



x4 E {23 + 1)1 , ^ 



2N 



S=-^N 



N-S 



exp 



2 M 



(3 J) {2S+lf 



We rewrite the last exponential by applying 



exp 



for p = 25* + 1, q = jPJ, and obtain 



Z = exp 



--a{N + l) 



J- 



dt 



dt exp 



1 9 

2q 



'^^2^a/{N + 1) 



X exp 



2a 



-\ 1^- I AT 



S=-^N V 2^ 



(59) 



(60) 



(61) 
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where we have introduced 



a = -l3J{N + l) 



(62) 



Now the summation over S can be carried out, leading to 



Z = exp 



-AN + l] 



dt 



^2TTa/{N + 1) 



exp 



2a 



d_ 

di 



expt cosh^t^ . (63) 



Upon using the t ^ —t symmetry of the integrand, introducing an additional derivative, and 
performing partial integration, we get 



Z = exp 



1 



a{N + l] 



dt 



2na/{N+l) 



1 f)"^ 

1 7V-I-1 



1 

'2^ 



{N + l)t' 



(64) 



With the equality 



a2 exp \-±{N+l)i 



p-na/iN + 1) 



this turns into 
Z = 2 exp 

= exp 



ra(A^ + 1) 



'-a{N + 1) 



d 



da 



2{N+1) 



dt 



d exp 



da 



^2Tia/{N + l) 



(65) 



27ra/(A^ + 1) 



exp 



— t^ 


cosh t^ 


2a 





N+l 



(66) 



°° (A^ + l)dri 
-^^2Ti/a{N + 1 



■■ T] taxih. arj exp 



-arj 



N+l 

cosh a?7 ) . (67) 



(Eq. (^) is obtained by computing 8/ da after substituting t t'^/a in eq. 

In the antiferromagnetic case, J < 0, the derivation is the same up to eq 
eq. (pup we then apply 



exp 



[ 1 21 


1 




1 2 • 




\/2'Kq . 


/ dt exp 

'—^00 


- t^ + ipt 
2q 



Instead of 



(6^ 



and the remainder of the derivation proceeds analogously to the ferromagnetic case, leading to 



-2 exp 



a{N + l] 



d_ 

da 



dt 



pT:a/{N + l] 



exp 



2a 



cost 



N+l 



(69) 



-a, t it. 



(where now a = —^[3J{N + 1) > 0). This expression is related to (|66D by a - 

For completeness we also mention the inclusion of an external field. If we add a coupling 
term —BS3 to the Hamiltonian, the factor (25" + 1) in eq. (|57| ) is replaced by 



^ ^ , sinh i/?S(2^ + 1) 



(70) 
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Proceeding along the lines of the derivation given above one obtains (with to 



z 



exp -\a{N +1) 1 




sinhio + 1 

d dt 



exp ——t^ cosh (t + to) 



) 



(71) 



^^0 J-^ ^27ra/{N + 1) 



in the ferromagnetic case, with a similar result for the antiferromagnet. 

References 

[1] I. Affleck, "Field Theory Methods and Quantum Critical Phenomena" (Les Houches lectures 
1988), in "Fields, strings and critical phenomena", eds. E. Brezin and J. Zinn- Justin (North- 
Holland, 1990) 

[2] J.H. Van Vleck, "The Theory of Electric and Magnetic Susceptibilities" (Clarendon Press, Oxford, 
1932) 

[3] P.W. Anderson, "Basic notions of condensed matter physics" (Benjamin, 1984) 
[4] S.K. Ma, "Statistical Mechanics" (World Scientific, Singapore, 1985) 

[5] E.P. Wigner, "Group theory and its application to the quantum mechanics of atomic spectra" 
(Academic Press, 1959), ch. 13 

[6] G. de B. Robinson, "Representation Theory of the Symmetric Group" (Edinburgh University 
Press, 1961) 

[7] H.E. Stanley, "Introduction to Phase Transitions and Critical Phenomena" (Oxford University 
Press, 1971), pp. 91-93 

[8] C. Kittel and H. Shore, Phys. Rep. 138 (1965) A1165 

Figure captions 

Figure 1: The integrand in eq. ( p9|) for a < 1 (solid line) and a > 1 (dashed line). 

Figure 2: Heat capacity per spin for the isotropic model (A^ = 2000) as a function of a = 
i/3J(iV + 1) (adapted from ref. §). The dotted line is | - ^(a - 1), see eq. (H), the dash- 
dotted line is eq. (|5|): 0.61 + 0.3672000 (a - 1). 
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